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ABSTRACT 

Using Suzaku and the Rossi X-ray Timing Explorer (RXTE), we have conducted a series of four simultane- 
ous observations of the galactic black hole candidate Cyg X- 1 in what were historically faint and spectrally hard 
"low states". Additionally, all of these observations occurred near superior conjunction with our line of sight 
to the X-ray source passing through the dense phases of the "focused wind" from the mass donating secondary. 
One of our observations was also simultaneous with observations by the Chandra-High Energy Transmission 
Grating (HETG). These latter spectra are crucial for revealing the ionized absorption due to the secondary's 
focused wind. Such absorption is present and must be accounted for in all four spectra. These simultaneous 
data give an unprecedented view of the 0.8-300 keV spectrum of Cyg X-1, and hence bear upon both corona 
and X-ray emitting jet models of black hole hard states. Three models fit the spectra well: coronae with thermal 
or mixed thermal/non-thermal electron populations, and jets. All three models require a soft component that 
we fit with a low temperature disk spectrum with an inner radius of only a few tens of GM/c^. All three models 
also agree that the known spectral break at 10 keV is not solely due to the presence of reflection, but each gives 
a different underlying explanation for the augmentation of this break. Thus whereas all three models require 
that there is a relativistically broadened Fe line, the strength and inner radius of such a line is dependent upon 
the specific model, thus making premature line-based estimates of the black hole spin in the Cyg X-1 system. 
We look at the relativistic line in detail, accounting for the narrow Fe emission and ionized absorption detected 
by HETG. Although the specific relativistic parameters of the line are continuum-dependent, none of the broad 
line fits allow for an inner disk radius that is > 40 GM/c^. 

Subject headings: accretion, accretion disks - black hole physics - radiation mechanisms: non-thermal - X- 
rays:binaries 



1. INTRODUCTION 

There is currently significant debate as to the physical 
mechanisms responsible for the continuum of X-ray spec- 
trally "hard states" of black hole candidates (BHC) accreting 
in binaries. This debate ranges from the broader issue of 
whether or not there is a significant contribution to the X-ray 
band from an outflow or jet, to more narrowly focused issues 
within given classes of models. For instance, the hard X-ray 
emission has traditionally been attribu ted to a Comptoniz- 
ing thermal corona (Eardle v. Lightm an & Shaiprd 1 19751; 
Shapiro, Liahtman & Eardlev Sunvaev & Triimpea 

' [l979i:iDove et al..,199 8. and references therein). These earlier 
works generally favor a scenario where the corona lies central 
to a truncated outer thin disk. However, if the corona is 
driven outwards by radiative pressure (e.g., Beloborodov 
1 19991) could it instead overlay the inner disk? Can the 
optically thick, geometrically thin dis k extend inward ne arly 
to the innermost stable circular orbit (iMiller et al.ll2006|) ? Is 
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this disk cold (peak temperatures of a few hundred eV), o r 
can it instead be hot (near a keV, i.e., I Wilms et al.l 120061) ? 
Is the hard state corona co mprised primarily of elec trons 
with a thermal population jPoutanen & Vumil 120091 ). or 
can it have a substa ntial contr ibution fro m a non-thermal 
electron population (Ibra gimov et al.ll2005h ? Does the bulk 
motio n of the flow play a role in Comptonizing the spec- 
trum (IShaposhnikov & TitarchukI 120061 : iLaurent & TitarchukI 
120071) ? Alternatively could the X-rays be comprised of a com- 
bination of direct synchrotr on and synchrotron self-Compton 
(SSC) emission from a jet (Mai'koff, Nowak & Wilriisl 120051: 
iMaitra e t al. 2009)? Contributing to the debate, however, is 
the fact that many of the above cited models, especially when 
considering solely Rossi X-ray Timing Explorer (RXTE) 
data in the 3-200 keV range (or even narrower energy 
ranges), describe the hard state spectra nearly equally well. 

To study these issues, over the past decade we have been 
using a series of pointed, approximately bi-weekly, RXTE 
observations of the BHC Cyg X-1 performed simultane- 
ously with 15 GHz radio observations by the Ryle telescope. 
Cyg X- 1 holds much promise for exploring the current range 
of questions listed above owing to its persistently bright X- 
ray flux (including both the hard and soft states, it varies be- 
tween 200-600mCrab in the 1.2-12 keV band covered by 
the i^TE - All Sky Monitor) a nd its correlated radio/X-ray 
spectra (see iWilms et alj|2006l and references therein). Ex- 
tended radio emiss ion even has been imaged in Cyg X-1 
(IStirling et al.ll200l . This campaign has already provided the 
spectr a for some of the Comptonization models (Wilms et al] 
120061) . and jet models jMarkoff. Nowak & Wilmsll2005i) dis- 
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cussed above. Additionally, these data have been used to 
to study the correlation of the X -ray and radio spectral 
properties on both long time scales ("Pottsch midt et al.ll2000t : 
[Ppttschmidt et al. 2003; Wilms et al. 2006) , and short time 
scales (iGleissner et all 120041 iWilmset alJ 12007. Boc k et all 
120 lOL Bock et al., in prep.). Thus they comprise a strong 
data set for addressing not only the details of the X-ray spec- 
trum, but also the connection to the hard state jets which are 
known to dominate the radio through near-infrared emission. 
Simultaneous radio/X-ray flaring has also been detected in 
Cyg X-1 during this extended campaign (£enderet al. 200^ 
I Wilms et alJ2067l) . lending support to the hypothesis of X-ray 
emission by the jet. 

RXTE spectral data alone, however, do not allow us to 
break the current existing "theoretical degeneracy" in the ori- 
gin of the X-ray spectrum. The statistically best fits are in 
fact obtained with purely empirical simple broken powerlaws 
with a break occurring between 9-12 keV, and an exponen- 
tial cutoff occurring at >20keV, to which a broad, » 6.4 keV 
gaussian fine is added jWilms et al.ll2006HNowak et alj2005h . 
The latter component is likely attributable to a rela tivistically 
broadened Fe Ka line ("Reynolds & NowalO 12003. and refer- 
ences therein); however, its para meters are d ependent upon 
the assumed continuum model (IWilms et alJ l2006l) . More 
physical ly motivated Comptonization and outflow-dominated 
models dMarkoff. Nowak & Wilmsll2005l: IWilms et al.ll2006l) 
can describe the same spectra almost as well. However they 
must introduce additional, albeit plausible, physical compo- 
nents (e.g., relativistic smearing) to recover the simple spec- 
tra of the broken powerlaw description. Thus, there is some 
amount of ambiguity when correlating detailed spectral fea- 
tures vs. continuum properties, e.g., disk reflection vs. coro- 
nal compactness/hardness, as the detailed features systemati- 
cally depend upon the underlying broad-band continuum. 

Despite these problems in finding a truly unique spectral 
model, when considering multiple observations taken over a 
wide range of luminosities and spectral hardnesses, spectral 
correlations arise that are robust and pers istent across a vari - 
ety of these theoretical characterizations ( IWilms et al.ll2006h . 
Using the broken powerlaw models as a simple description 
of the X-ray spectra, we have found that when the 2-lOkeV 
photon index Fi < 2.2 there is a positive correlation between 
X-ray and radio flux, whereas for Fi > 2.2 there is a negative 
correlation between the X-ray and radio flux. We therefore 
use the value of Fi « 2.2 as the canonical division between 
the spectrally "hard state" and the spectrally "soft state" (see 
iRemillard & McClintock 2006). Additionally, as the spectra 
become harder, exponential cutoffs tend to beco me less sig- 
nificant (see also lMotta. Belloni & Homanll2009l) . 

Our previous RXTE spectral studies of Cyg X-1 have been 
limited in two respects: the low spectral resolution of RXTE 
(E I A/i ss 6 at 6 keV), and the inability to measure spectra at 
^ 3 keV. In this work, we turn to a set of four Suzaku ob- 
servations that we performed simultaneously with our RXTE 
campaign to enhance our Cyg X-1 studies in several crucial 
ways. First, Suzaku has large effective area at soft X-ray en- 
ergies. In this work, we consider spectra down to 0.8 keV, 
which allows for the possibility of measurin g the "seed pho- 
ton spectrum" in Comptonization models 05.11 ). or judging 
the relative contribution of synchrotron versus the disk radia- 
tion in jet models 05.2l i. Suzaku also has excellent resolution 
in the Fe Ka line region (E/AE w 50), which allows us to 
s eparate na rrow from relativistically broadened line features 
(£3 |53]l. Third, Suzaku measures the Cyg X-1 hard X- 



TABLE 1 
Log of Cyg X-1 Observations 



Date Spacecraft/ObsID Instrament Exposure 

(yyyy-mm-dd) (ksec) 



2006-10-30 


SuzaJcu/40 10590 10 


XIS 0-3 


35.0" 


32.1 






HXD-PIN 


27.7 


24.9 






HXD-GSO 


27.7 


25.8 




RXTB/801 10-01-13 


PCA 


8.0'' 


8.0 






HEXTE-A 


3.1 


3.1 






HEXTE-B 


2.5 


2.5 


2007-04-30 


Suzafcu/402072010 


XIS 0,1,3 


34.0" 


22.4 






HXD-PIN 


40.2 


27.6 






HXD-GSO 


40.2 


27.1 




RXTE/92090-01-16 


PCA 


14.6'' 


1.5 






HEXTE-A 


4.6'^-'' 


0.3 






HEXTE-B 


3.0 


0.4 


2007-05-17 


SuzaJcu/402072020 


XIS 0,1,3 


22.3" 


12.3 






HXD-PIN 


32.6 


10.9 






HXD-GSO 


32.6 


10.9 




RXTE/92090-01-17 


PCA 


2.0'' 


0.0 






HEXTE-A 


1.9'^'' 


0.0 






HEXTE-B 


0.9 


0.0 


2008-04-19 


Suzafcu/403065010 


XIS 1,3 


16.9" 


7.2 






XISO 


34.0= 


0.0 






HXD-PIN 


29.0 


10.1 






HXD-GSO 


29.0 


10.1 




RXTE/93 120-0 1-01 


PCA 


21.5'' 


7.2 






HEXTE-A 


17.9"^ 


7.1 






HEXTE-B 


10.8 


4.6 




CJiandra/8525 


HETG 


29.4'' 


11. IS 



Note. — Exposure times ai'e after initial good time filtering (left), 
and after color/intensity time filtering (right). 
" Summed exposure times for all listed XIS detectors. 
'' Summed exposure intervals, not weighted by the fraction of operating 
Proportional Counter Units (PCU). 

HEXTE-A cluster in fixed position (no rocking). 

Evidence for HEXTE-A cluster exposure time anomaly. 

XIS run in continuous readout mode. 
^ Summed exposure times for all l"' order spectra. 
^ Summed exposure times for 1^' order HEG spectra. 

ray spectrum up to w 300 keV 04.1l i. providing further con- 
straints on Comptonization and jet models. 

The outline of this paper is as follows. In ^ we describe 
our full set of observations and our data reduction procedures. 
Due to maintenance and upgrade of the Ryle radio telescope 
during the construction of the Arcminute Microkelvin Imager 
(AMI), no simultaneous radio measurements are available for 
these SuzakulRXTE observations. Inste ad, w e use previous 
observations to estimate the radio fluxes 02.4| i. For one of our 
observations simultaneous Chandra -HETG data are available 
02.3I ). These Chandra data become crucial in all of our anal- 
yses as they help elucidate the spectral variability associated 
with the observed lightcurve behavior, as discussed in ^ We 
present simple phenomenological descriptions of the spectra 
i n m including a description of the composite line profile 
04.3l l. Comptonization and jet models are presented in ^ 
al ong with further discussions of the irnplied relativistic lines 
05. 3b . We summarize our findings in ^ 

2. OBSERVATIONS AND DATA ANALYSIS 

Prior to 2009 April, there have been five Suzaku observa- 
tions of Cyg X-1. The four most recent of these observations 
are discussed in this work. The fi rst SuzafcuCyg X-1 ob - 
servation occurred in 2005 October (iMakishima et al.|[2008h . 
During this observation Cyg X-1 was in a relatively bright 
hard state (the RXTE-ASM count rate was « 30cps), the 
spacecraft aimpoint was placed on the Suzatu -X-ray Imag- 
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ing Spectrometer (XIS) detectors, and Suzaku was run in a 
data-mode w ith 1 s integrations per CCD exposure frame. As 
discussed bv'M akishima et all ( 12008 , see also the Appendix), 
this long exposure led to both telemetry dropouts and severe 
photon pileup on the detectors. As these issues require a more 
complex analysis (Makis hima et al. 2008 ), we do not consider 
these data further. 

The four Suzaku observations discussed here occurred dur- 
ing times when the RXTE-ASM flux ranged from 12-23 cps, 
the spacecraft pointing was set to the Hard X-ray Detector 
(HXD; this slightly reduces the flux on the XIS detectors), 
and the CCD exposure frame integration times were set to 
« 0.5 s. Thus, our observations do not suffer from telemetry 
dropouts and they are less severely affected by photon pileup. 
All four of these observations occurred simultaneously with 
RXTE observations. The last of these observations also oc- 
curred simultaneously with observations by every other X- 
ray/soft gamma-ray instrument flying at that time (Hanke et 
al., in prep.). Here we only consider the simultaneous Chan- 
dra-High Energy Transmission Grating. An observing log is 
presented in Table [T] 

2.1. Suzaku Analysis 

The Suzaku data were reduced with tools from the HEA- 
SOFT v6.8 package and calibration files dat ed 2009 Sept em- 
ber 25. The instruments on Suzaku ( Mitsud a et al.ll2007h are 
the X-ray Imaging Spectrometer (XIS: iKovama et al.l 120071) 
CCD detector covering th e w 0.3-lOkeV band , and the Hard 
X-ray Detector (HXD; Takahashi et al.l 12007 *) comprised of 
the PIN diode detector (PIN) covering the « 10-70keV band 
and the gadolinium silicate crystal detector (GSO) covering 
the « 60-600 keV band. The XIS has four separate detec- 
tors, XIS 0-3, with XIS 1 being a backside illuminated CCD. 
XIS 2 was lost due to a micrometeor hit in late 2006, and thus 
was available only for the first observation (see Table[T]). For 
our fourth observation, XIS was run in continuous readout 
mode, which is not yet fully calibrated, so we do not include 
these data in this work. 

In preparing the XIS spectra, we first corrected each detec- 
tor for Charge Transfer Inefficiency (CTI) using the xispi 
tool, and then reprocessed the data with xselect using the 
standard xisrepro selection criteria. Due to thermal flex- 
ing of the spacecraft, the attitude of the Suzaku spacecraft 
exhibits variability over the course of the observations and 
therefore the image of the source is not at a fixed position on 
the CCD. Standard processing reduces this variability and im- 
proves the PSF image (Uchiyama et al. 2008); however, with 
the standard tools it is not possible yet to correct fully the 
blurring caused by the varying attitude. As described in the 
Appendix, for bright sources a better reconstruction of the 
attitude solution, and thus narrower PSF images, can be ob- 
tained using the aeattcor.sl software. We produced such 
an improved image and then used the pile_est imate . si 
S-Lang script described in the Appendix to estimate the de- 
gree of pileup in the observations. For the spectra described 
in this work, using the time filter criteria described in ^ the 
center of the PSF images could be affected by as much as a 
35% pileup fraction. We therefore extracted annular regions 
wherein we excised the w 20" radius central region. The ex- 
cised data accounted for approximately 1/3 of the detected 
events. The outer radii of our annular extraction regions were 
limited by the 1/4 sub-array used in our CCD readout mode, 
and thus were w 2'. We estimate that the extracted events had 
< 4% mean residual pileup fraction. All regions of the CCD 



are dominated by source counts, therefore we did not extract 
nor use any background spectra for the XIS observations. 

Events in the XIS detectors were read out from either 3x3 
or 2 X 2 pixel islands. The XIS 1 data were always in 3 x 3- 
mode, whereas the other detectors had mixtures of 3 x 3- 
and 2 x 2-mode. We created individual spectra and response 
files for each detector and data mode combination. Spec- 
tra were created for the time intervals described in ^ with 
the resulting exposure times also being listed in Table [T] Re- 
sponse matrices and effective area files were created with the 
xisrmf gen and xissimarf gen tools, respectively. 

Although for each observation the individual XIS spectra 
were fit separately, they were jointly grouped on a common 
grid such that they had a minimum combined signal-to-noise 
ratio of 8 in each energy bin (i.e., 64 total counts in each bin) 
and that the minimum number of channels per energy bin was 
at least the half width half maximum of the spectral resolu- 
tiorQ To avoid regions of poorly understood response, we 
only considered spectral energy ranges 0.8-1. 72 keV, 1.88- 
2.19keV, and 2.37-7.5 ke\0 

For the PIN, we extracted spectra from the "cleaned" 
event files in the hxd/event_cl directories. The ap- 
propriate response and background files were downloaded 
from the High Energy Astrophysics Science Archive Re- 
search Center (HEASARC), specifically those from the 
pinxb_ver2 . 0_tuned directory. Good time intervals 
(GTI) were merged from our time selections (® and the 
combination of the GTI intervals from the source and back- 
ground event files. These intersected intervals were then 
used to extract the PIN source and background spectra. The 
source spectra exposure times were then corrected with the 
hxdtcorr tool. The PIN spectra were grouped to have a 
signal-to-noise ratio > 10 in each energy bin, and we consid- 
ered spectra between 12-70keV. 

The GSO spectra were created starting with the "unfil- 
tered" event files. These data were first reprocessed with 
the hxdtime, hxdpi, and hxdgrade tools, following 
the "Suzaku ABC Guide". The data were then filtered in 
xselect with the standard criteria from the HEASARC pro- 
vided gso_mkf .sel script. The background was down- 
loaded from the gsonxb_ver2 . directory at HEASARC. 
GTI from the event file, the background file, and the time in- 
tervals were merged, and spectra were extracted from these 
times. Response files were then taken from the CALDB 
database, and exposure times were adjusted to agree with the 
spectra. The grouping of the GSO spectra is essentially fixed 
by the grouping of the background file, thus no rebinning was 
performed on these spectra. It was determined that the back- 
ground becomes prohibitive above 300 keV; therefore, we re- 
strict the GSO spectra to the 60-300 keV range. 

2.2. RXTE Analysis 

' The half width half maximum (HWHM) was determined by using the 
spectral responses to create fake spectra of delta function lines at discrete 
energies. These fake spectra were then fit without response matrices to de- 
termine HWHM. In practice, this meant that the original spectra with 4096 
channels each were grouped, using the ISIS group function, by a minimum 
of 6, 8, 12, 14, 16, 18, 20, 22 channels starting at 0.5, 1, 2, 3, 4, 5, 6, and 
7 keV, respectively. For the most part, the spectral binnings were dominated 
by the HWHM criteria rather than the signal-to-noise criterion. 

^ Outside of these energy ranges there are large disagreements among the 
fit residuals for the individual detectors and data modes. The excised regions 
between 1.72-2.37 keV correspond to poorly calibrated Si and Ir features re- 
lated to the detectors and mirrors. 
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The RXTE data were prepared with tools from the HEA- 
SOFT v6.8 package and the most current caUbration files as 
of 2010 January 10. We used standard filterin g criteria for 
data f rom the Proportional Counter Array (PCA: Ijahoda et all 
[1996'). Specifically, we excluded data from within 30 minutes 
of South Atlantic Anomaly (SAA) passage, from whenever 
the target elevation above the limb of the earth was less than 
10°, and from whenever the electron ratio (a measure of the 
charged particle background) was greater than 0. 15. We used 
the background models appropriate to bright data. 

We applied 0.5% systematic errors to all PCA channels, 
added in quadrature to the errors calculated from the data 
count rate. For all PCA fits we grouped the data starting at 
> 3 keV with the criteria that the signal-to-noise (after back- 
ground subtraction, but excluding systematic errors) in each 
bin had to be > 4.5. We restricted the noticed energy range to 
3-22 keV. 

We extracted data f rom the High Energy X-ray Timing Ex- 
periment (HEXTE; Rothschild etalJl998|) using the same cri- 
teria as for the PCA. HEXTE is comprised of two clusters, A 
and B. Prior to 2007, both clusters were in a rocking mode, 
with only one detector viewing the source at a given time, 
while the other detector conducted off-source background 
measurements. For later dates (e.g., our three most recent ob- 
servations), HEXTE A has been in a fixed on-source position, 
which required us to estimate cluster A backgrounds using 
the hextebackest tool. Furthermore, during our third ob- 
servation there was clear evidence for a HEXTE A exposure 
anomaljQ We found evidence for such an anomaly (albeit 
less severe) in our second observation as well. We corrected 
the HEXTE A exposure time by iteratively fitting an exponen- 
tially cutoff powerlaw simultaneously to both cluster A and B 
data and varying the exposure of cluster A until a minimum 
was achieved using the same normalization for both detectors. 

We binned the HEXTE spectra to a common grid such that 
they had a minimum combined signal-to-noise of 8 in each en- 
ergy channel (even though both clusters were fit individually). 
We further restricted the noticed energy range to 1 8-200 keV. 
Additionally, when fitting the X-ray spectra we allowed the 
normalization of the HEXTE backgrounds to vary. (The best 
fit normalization constants were typically within < 10% of 
unity.) 

For both PCA and HEXTE spectra we further restricted the 
considered time intervals to those that were strictly simultane- 
ous with the Suzaku spectra. That is, the RXTE time intervals 
form a subset of the Suzaku time intervals. 

2.3. Chandra-HETG Analysis 

The High Energ y Transmission Grating (HETG; 
iCanizares et al.l I2005h was inserted for our observation 
of Cyg X-1, with the data readout mode being Timed 
Exposure-Graded. The HETG is comprised of the High 
Energy Gratings (HEG), with coverage from « 0.7-8 keV, 
and the Medium Energy Gratings (MEG), with coverage 
from « 0.4-8 keV. To minimize pileup in the gratings spectra, 
a 1/2 subarray was applied to the CCDs. Additionally, the 
observatory aimpoint was placed closer to the CCD readout. 
This configuration reduces the frame time to 1.741 sec, 
without any loss of the dispersed spectrum. 

We used CIAO v4 . 2 and CALDB v4 . 2 . to extract 
the data and create the spectral response files. The loca- 

^ http://gsfc.nasa. gov/ doc s/xte/ what snew/ 

newsarchive_2007 .html/ 




Time (sec) 

Fig. 1. — Suzatu-XIS 1 lightcurves (16 sec integration bins) for observa- 
tion 4 in the 0.5-1.5 lieV band (filled diamonds, bottom curve) and the 3- 
9keV band (hollow diamonds, top curve, shifted upward by 50cps). 

tion of the center of the 0* order image was determined us- 
ing the findzo.sl routing which provides w 0.1 pixel 
(« 0.001 A, for MEG) accuracy. The data were reprocessed 
with pixel randomization turned off, but pha randomization 
left on. We applied the standard grade and bad pixel file fil- 
ters, but we did not destreak the data. 

Although the instrumental set up was designed to mini- 
mize pileup, it is still present in both the MEG and, to a 
lesser extent, the HEG spectra. We incorporate the effects 
of pile up in our sp ectral fits using the simple_gpile2 
model (Hanke et al. 2009) based upon the model originally 
described by Nowak et al. (2008). For the small amount of 
pileup present in these data, this model is extremely success- 
ful in describing its effect on the spectra. 

As for the PCA and HEXTE spectra, we restricted the con- 
sidered time intervals for the Chandra spectra to those that 
were strictly simultaneous with the Suzaku spectra. 

2.4. Estimated Radio fluxes 

Simultaneous Ryle radio fluxes are not available for the X- 
ray observations discussed in this work. Owing to the known 
stringent correlations between the X-ray and radio flux, 
specifically the correlations between the RXTE -ASM and 
Ryle fluxes (see lGleissner et alJ2004tlNowak et alJ2005 h. we 
use an estimated radio flux when applying the jet model to 
the spectra 05.2l i. From our prior observations we average 
all radio data that occurred at times with ASM daily aver- 
ages that were within ilcr of the ASM daily average corre- 
sponding to these new observations. The average and stan- 
dard deviation of these radio data were used in the spectral 
fits, which yielded radio fluxes of 12.6 ± 3.3 mJy (2006-10- 
30), 10.8±3.9mJy (2007-04-30), 7.9±3.0mJy (2007-05- 
17), and 10.2±3.0mJy (2008-04-19). Each average was 
comprised of a minimum of ten measurements, and the ex- 
treme radio flux values within each sample never differed 
from the average by more than a factor of two. 

2.5. Data Plots 

http : / / space . mit . edu/ASC/ analysis/ f indzo/ 



Observations of the Cygnus X-1 Hard State 
in 




0.4 0.6 

(0.5-1.5 keV)/(1.5-3keV) 



0.2 0.4 0.6 0.8 

(0.5-1.5 keV)/(1.5-3keV) 



Obs. 3 




in 



On 
I 



> 

m 

I 



in 
o 



Obs. 4 




0.2 



0.4 



0.6 



0.8 



(0.5-1.5 keV)/(1.5-3keV) 



(0.5-1.5 keV)/(1.5-3keV) 



Fig. 2. — Suzatu-XIS 1 color-color diagrams, made with 16s bins. The colors are created from the count rate in three bands: 0.5-1.5 keV, 1.5-3keV, and 
3-9 keV. The lines are the best fit dust halo models, presuming that the source is obscured by both an interstellar neutral column as well as by a neutral column 
local to the source, with the latter ranging in value from 0-2 X 10^' cm"^. Filled symbols indicate the periods used for the spectral analyses described in this 
work. (See text.) 



Throughout this work we will present spectra that are "un- 
folded" fro m the detector respo nse in a model-independent 
manner (see 'N owak et al.ll2005l) . Specifically, we define the 
unfolded photon flux, Faafo\A{h), from a pulse height (PHA) 
bin, h, as 



i^unfold 



{h): 



C(h)-B(h) 



J R(h,E)A(E)dEAt ' 



(1) 



where C{h) - B(h) are the background subtracted counts, 
R(h,E) and A{E) are the detector response matrix and effec- 
tive area, and At is the exposure time. For multiple data 
sets coadded together for the plots (e.g., the Suzaku-XlS or 
the RXTE-HEXTE data), the numerator and denominator are 
each replaced with the sum of the same quantities from the 
detectors. Throughout we only sum data for the plots, not for 
the fits. 

Model spectra are "unfolded" in a similar manner. That is, 
predicted background subtracted model counts are first calcu- 
lated, and then unfolded as above. Thus, the model spectra 



exhibit the smearing induced by the detector response matrix. 
The one exception to th is m odel component scheme is for the 
jet results presented in ^5.21 There some of the model compo- 
nents are shown unsmeared by the detector response, as they 
are derived from internal calculations of the jet code. 

3. LIGHTCURVES 

By design, our fourth observation occurred near binary or- 
bital phase (superior conjunction of the black hole), wherein 
we expected to be viewing Cygnus X-1 through the wind of 
the secondary. By happenstance, all of our remaining Suzaku 
observations also occurred near superior conjunction, i.e., ob- 
servations 1-4 cover orbital phases 0.2-0.3, 0.8-0.0, 0.85- 
0.0, and 0.0-0.14, respectively. Although dipping events are 
known to be present during all orbital phases, during upper 
conjunction dipping is prevalent due to absor ption by clumps 
in the wind (iBalucinska-Church et alj|2000l) . This behavior 
is clearly seen in the Suzaku lightcurves presented in Fig. [T] 
The soft X-ray lightcurve shows a deep, prolonged « 12ksec 
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dip. A less pronounce d, but still significa nt dip is seen in the 
hard X-ray lightcurve dHanke et al.ll2008h . Numerous dips of 
varying levels are also seen on even shorter time scales. In 
order to extract the unabsorbed spectrum of the source, it is 
therefore necessary to screen the data for dipping and excise 
the dipping intervals from further analysis. 

In addition to dipping, spectral modeling also has to 
take into account the strong dust scattering halo in front 
of Cyg X-1 jPr edehl & Schmitt 1995). As discussed by 
IXu. McCray & Kelleyi (il986t and references therein), a scat- 
tering halo is produced by the scattering of radiation from an 
X-ray source by a foreground dust cloud. Typical halo sizes 
are on the order of arcminutes, i.e., comparable to the size 
of the Suzaku PSF. The energy dependence of the scattering 
cross section leads to a halo angular size that is oc and 
a halo spectrum that is oc E~^. Furthermore, as the scattered 
photons travel along a greater path length to reach the ob- 
server, they are delayed with respect to the directly observed 
photons by a factor that is proportional to the square of the 
angular radius from which the scattered photons are observed. 
The expected time delays for nearby galactic sources typically 
range from thousands to tens of thousands of seconds. 

Thus the Suzaku spectrum and lightcurve of Cyg X-1 are 
comprised of two components: a directly observed compo- 
nent, subject to the local dipping events, and a time-delayed 
and time-averaged (due to integration over different angular 
radii) softer spectrum from the scattering halo. Most spec- 
tral analyses in the past have ignored the effects of scattering 
since to first order for optically thin scattering in a homoge- 
neous cloud, whatever radiation that is scattered out of the 
line of sight is scattered back in from larger (often spatially 
unresolved) radii. Thus, for spatially unresolved (on the size 
scale of the halo) and steady sources, dust can be ignored in 
simple analyses. 

If one assumes that the time-delayed and time-averaged 
spectrum is represented by the same model as the direct spec- 
trum, and if one further assumes that interstellar absorption 
is predominantly in the foreground of the dust halo, within 
ISIS one can write the spectral model for the spectra during 
the dips as: 

TBnew( 1) * ( 1 - dust s cat * ( 1 - TBnew(2)))* 

(contirLuum model) . (2) 

Here and throughout we shall use an updated v e rsiorFI of 
the absorption model of 'Wilms. Allen & McCravl (|2000j to 
describe both the interstellar (model instance 1) and local 
(model instance 2) absorption. The dustscat model is a 
version of the XSPEC dust model that removes the opti- 
cally thin assumption of the latter model, but otherwise ig- 
nores multiple scatterings (F. Baganoff, priv. comm.). Note 
that when local absorption is absent (TBnew (2 ) 1), the 
dust scattering term drops out of the model expression. 

Not all of the assumptions encompassed within the above 
model expression are necessarily realistic. Specifically, the 
assumptions that all interstellar absorption is foreground to 
the dust and that the delayed source spectrum is identical to 
direct spectrum are undoubtedly wrong in fine detail. How- 
ever, we find that the above expression provides a reasonable 
description of the color-color diagrams in Fig.|2l Specifically, 
we create a spectrum from bright phases of the lightcurve that 

^ http : //pulsar . sternwarte . uni-erlangen . de/wilms/ 
research/tbabs 



TABLE 2 

Observational Fluxes in keV Energy Bands 



Date 


0.5-2 


2-10 


10-100 


100-300 


Lfiol" 


(yyyy-mm-dd) 




(10"** 


erg cm"^ s" 




(Lfidd) 


2006-10-30 


0.16 


0.78 


2.76 


1.59 


0.032 


2007-04-30 


0.07 


0.55 


2.07 


1.32 


0.026 


2007-05-17 


0.07 


0.38 


1.33 


0.87 


0.017 


2008-04-19 


0.14 


0.72 


2.62 


1.44 


0.035 



Note. — Fluxes are absorbed values, are normalized to the 
PCA spectral fit, and correspond to the brightest/least absorbed 
periods of the lightcurves. 



" Unabsorbed, isotropic luminosity in 0.01-800 keV band, ex- 
pressed as a fraction of Eddington lu mino sity, as determined from 
the thermal Comptonization fits of j|5.1l and assuming a distance 
of 2.3 kpc and a black hole mass of IOMq . 

inhabit a locus of points in the upper right hand corner of the 
color-color diagrams (i.e., the filled symbols of Fig. |2|l. We 
fit a simple spectral model to these selected data: an absorbed 
disk plus powerlaw plus broad and narrow gaussian lines near 
6.4 keV. We then apply this model in eq. (|2]i, and create color- 
color curves by varying the TBnew ( 2 ) component from 0- 
2 X lO^^cm"^. These curves are then fit to the color-color 
diagrams using the optical depth of the dust scattering halo as 
the single free parameter Our best fits yield r = 0.24-0.34, 
and the results are shown in Fig.|2] 

Overall, these curves describe the behavior of the color- 
color diagrams reasonably well. The dust halo represents 
« 20%-30% (i.e., on the order of the scattering optical depth) 
of the soft X-ray (< 2 keV) flux being "uncovered" during 
the dips. Given the much narrower Chandra PSF, one would 
expect comparable Chandra color-color diagrams to show a 
far smaller uncovered fraction if the dust halo interpretation 
is correct. This is indeed found to be t rue, with Chandra 
showing only an « 2% uncovered fraction (H anke et al.ll2008l 
Hanke et al., in prep.). However, it is possible that the very 
low uncovered fraction required by the Chandra data might, 
in fact, be partly intrinsic to the local dipping at the source 
(Hanke et al. in prep.). Note that the theoretical curves show 
a greater degree of curvature than the observed, nearly linear, 
evolution toward the lower left hand corner of the diagram. 
We hypothesize that this discrepancy is due to our assumption 
in eq. (|2]) of purely neutral absorption in the dips. A concur- 
rent increase in the optical depth of ionized absorption could 
serve to alter the theoretical curves in the appropriate manner 
We return to this concept of an ionized absorber component 
inS2] 

4. SIMPLE SPECTRAL MODELS 

The lightcurves and color-color diagrams described above 
were used to create spectra for the time periods corresponding 
to the bright phases of the lightcurve represented in the upper 
right corner of the color-color diagrams (Fig. |2]l. The selec- 
tions that we made are highlighted as filled symbols within 
these figures. The resulting exposure times are given in Ta- 
ble[T] With the exception of the RXTE-only broken powerlaw 
spectral fits described immediately below, throughout the rest 
of this work we shall only describe spectra from these bright 
phase, minimaUy locally absorbed time periods. 

4.L Broken and Cutoff Powerlaw Descriptions 

As discussed by ' Wilms et all (120061) , nearly afl RXTE- 
PCA and HEXTE spectra of Cyg X-1 can be described 
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TABLE 3 

CONSTANT*! I GHE CUT *POWERLAW PARAMETERS FOR HARD X-RAY SPECTRA FITS 



Date 
(yyyy-mm-dd) 


r 


ApL 


(keV) 


Efold 

(keV) 


Cgso 


Chexte a 


Chexte b 


nnorm 

''hexte a 


nnorm 

"hexte B 


X^/DoF 


2006- 10-30 

2007- 04-30 

2007- 05-17 

2008- 04-19 


J 4J+0.02 
^■^^-0.01 

1 47+0.02 

1 iq+S:8f 

'■-^ -0.01 


1.07+0°' 

7S+^:8^ 
"•"-0.03 

60+'"''* 
-Q.Q3 

"•"-0.03 


19.6+3-0 
24.6!|:» 

20A+\l 


183^^ 
218+^ 
252+2' 
164+° 


1 007+O-008 

'-""'-0.008 

95+0 01 

"-"-0.01 

95+0 02 
9135.01 

"-^'-0.01 


n 077+0.004 
1 n7+0-02 

'-"'-0.03 

n o^7+0.003 
O-'*^ /_o (,03 


n ss^+0.005 
1 07+0-02 

'-"'-0.01 

n S/19+0.004 
"-'^4^-0.004 


1 12+001 

1 03+^:'l'o 
'-"-"-0.06 

1 065+^-008 

1 ."u J_o 008 


1 02+002 
-8!) J 

n Qn+0.06 
"-^^-0.06 

99+0.01 

"-^^-0.01 


785.4/449 
966.2/315 
301.6/123 
648.8/447 



Note. — Spectra ai'e normalized to HXD-PIN data. E^vit, £fold are the cutoff and folding energy of the highecut model. Api is the 
power law normalization in photons keV"' cm"2 g-l at j j^gV. C are the fit constants for the detectors other than the HXD-PIN. The normalizations 
of the HEXTE backgrounds were also allowed to be free parameters, and are given by the B"™™ values. 



by an extremely simple phenomenological model: an ab- 
sorbed, exponentially cutoff, broken powerlaw plus a broad, 
gaussian line. The photon indices of the broken pow- 
erlaw, which we label as Fi and r2 for the soft (« 3- 
lOkeV) and hard {<L lOkeV) X-rays, show a correlation 
in the sense that the amplitude of the break between the 
two, Fi — increases with higher (i.e., softer) values of 
Fi. This F-AF correlation is a very phenomenological 
description of the spectral correlation that elsewhere has 
been descri bed as a "T-U/2tt" (hardness /reflection fraction) 
correlation (IZdziarski. Lubih ski & SmithI [T999i). Fo llowing 
iRemillard & McClintock (2Q06h. IWilms et all ( 120061) classi- 
fied spectra with Fi < 2.2 as "hard state" spectra. 

In order to compare the spectra discussed here to the work 
of IWilms et alj (2006), we take the 9 RXTE spectrtS and fit 
them wi th exponenti a lly cu toff broken powerlaws. As in the 
work of IWilms et alJ (120061) . these simple phenomenological 
models are excellent descriptions of these data. This is even 
true for the spectra extracted from times of moderate and 
deep dips. These latter spectra have slightly reduced values 
of Fi and increase d values of the fitted neutral column, A^h 
( IWilms et al.ll2006l show that orbital phase-dependence of the 
column is discernible in our RXTE observations of Cyg X-1), 
but otherwise appear similar to non-dip phase spectra. This is 
not entirely surprising given that RXTE is primarily sensitive 
at energies > 3 keV, and therefore does not measure the spec- 
tral regime most strongly affected by the dips. 

We show the resulting p hoton indice s of these spectra 
along side the results from IWilms et alJ (^06) in Fig. [3] 
Our prior study covered seven years and 202 RXTE spec- 
tra. These newer observations overlap the historically hardest 
states measured in that campaign. Specifically, for these new 
observations Fi ranged from 1.63-1.71, whi ch is to be com - 
pared to the Fi = 1 .65 minimum found in Wi lms et alj (120061) . 
In some cases, these observations also correspond to histori- 
cally faint hard states of Cyg X-1. Defining the "hard state" as 
Fi < 2.2, the lowest 2-100 keV flux in the hard state reported 
by IWilms et"an C2006) was 2.46 x lO'^^ergcm'^s"', while the 
highest 2-100 keV hard state flux reported in that work was 
4.50 X lO-'^ergcm-^s"'. As shown in TablelH the 2-lOOkeV 
fluxes spanned by the bright phase, least absorbed observa- 
tions range from 1.71-3.54x lO'^^ergcm'^s"'. 

That these spectra are among the faintest and hardest 
ever observed for Cyg X-1 is important to note. In some 
models of the black hole hard state, the hardest spectra 

^ These 9 spectra consist of three bright phase spectra, and six spectra 
taken from deeper dipping periods. Observation 1 shows no strong dips, and 
observation 3 does not have RXTE data strictly simultaneous with the bright 
phases of the Suzaku spectra. 




1.5 2 2.5 3 



r, 

Fig. 3. — Correlation of the soft X-ray (Fi) and hard X-ray (r2) photon 
indices from exponentially cutoff, broken powerlaw fits to the RXTECyg X- 
1 data presented by Wilms et al. 1 2006). Triangles represent the hard states 
(Fi < 2.2), and circles represent soft states (F2 > 2.2). The solid diamonds 
represent exponentially cutoff, broken powerlaw fits to the RXTE data dis- 
cussed in this work. There are spectra from four observations that have been 
further subdivided into up into times covering different portions of the color- 
color diagrams of Fig.|2] for a total of 9 RXTE spectra. 

are expected to have the smallest observed reflection frac- 
tion and the disk inner radius should be at a maximum 
dZdziarski. Lubinski & Sm ith 1999). However, it is also im- 
portant to note that other black hole systems can show sub- 
stantially fainter spectra relative to Eddington luminosity. 
Furthermore, as we see here and as has been noted elsewhere, 
the range of hard state fluxes observed in Cyg X-1 spans only 
a factor of two, and it is unlikely that the bolometric flux spans 
a much greater factor than this even when including observa- 
tions from the soft state (Wilms et al. 2006). 

As a further comparison to the previously observed 
hard states, we consider the value of the fi tted folding 
energy . For the hard states discussed by IWilms et alj 
(I2OO6I) . the folding energy ranges from ss 125-255 keV. 
Other BHC, e.g., GX 339-4, have shown a wi der range of 
hard state folding energie s, i.e. , «50-300keV (Wilms et alj 
1999t Nowak. Wilms & Dov3 l2002t .Nowaketal. 2005t 
Motta. Belloni & Homan 2009). In our prior studies we did 
not consider spectra above 125keV. Here, however, with the 
inclusion of the Suzaku-GSO spectra we can now consider 
energy ranges up to 300 keV. For all four of our observations, 
an exponential rollover is clearly detected. An example is 
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Energy (keV) 

Fig. 4. — Suzaiu-HXDPIN (brown squares) and GSO (magenta triangles) 
spectra, and RXTiJ-HEXTE (green diamonds) spectra from the least ab- 
sorbed periods of observation 4 fit with an exponentially cutoff powerlaw. 

given in Fig. 21 which also gives a general indication of the 
consistency of the cross-coiTelation among the hard X-ray de- 
tectors. The values of the folding energies, along with cross- 
normalization constants, are given in Table [3] The folding 
energies range from 164-252 keV, nearly the full span encom- 
passed by our prior hard state observati ons. This is not neces- 
sarily surprising. Whereas Wil ms et"an (f2006) had observed a 
general trend for the cutoff to increase to higher energies with 
harder spectra, there was a large degree of scatter about this 
trend, consistent with these four observations. 

4.2. Ionized Absorption Models 

Previous Chandra-HETG observations of Cyg X-1 have re- 
veal ed evidence of the ionized, focused w ind from the donor 
star (iMiller et al .1120051; iHanke et al.ll2009h . The effects of this 
wind upon the spectrum must be accounted for in any de- 
tailed analysis of Cyg X-1. While ionized line absorption is 
very pronounced during the deepest part of the absorption dips 
(Fig.[T]l, we stress that such absorption is also significant dur- 
ing the bright phase (non-dip) spectra discussed in this work. 

We begin by considering the joint 0.8-7.5 keV Suzaku 
and Chandra-RETG spectra from the 2008 April 19 ob- 
servations. As the MEG spectra are more heavily af- 
fected by photon pile-up, we fitted only the HEG spec- 
tra. We used the same absorbed disk plus powerlaw 
(TBnew* (diskbb+powerlaw) ) continuum model to fit 
both the Suzaku and HETG spectra. We also added to 
both spectral models a broad Fe Ka line (modeled with 
diskline ; iFabian et al.lll989h and a narrow Fe Ka line 
modeled with a gaussian, and a 6.7 keV (Fe XXV) absorption 
line. Following Han ke et al.. (^009), we describe the line ab- 
sorption as a series of gaussian absorption lines representing 
a variety of ionized species — predominantly H- and He-like 
lines from elements ranging from O to Fe, and L-shell transi- 
tions of Fe. This initial line list contained nearly 100 line tran- 
sitions parameterized by wavelength (which we constrained to 
fall within 1500 kms"' of their rest wavelengths), full width 
half maximum (constrained to be between 5^0 mA), and line 
equivalent width. Based upon fits to the joint Suzaku-HETG 
spectra, this initial list was reduced to the 55 significant lines 
present in the spectra. 



at 




Energy (keV) 

Fig. 5. — Suzaku (blue diamonds) and CAandra-HETG (grey histogram; 
iP'order HEG only) data from observation 4. The model consists of ab- 
sorbed disk and powerlaw components, broad and narrow gaussian lines in 
the Fe Ka region, and absorption lines from a variety of highly ionized 
species, including from Fe XXV at 6.7 keV. Additionally, we modify the 
HETG spectra by dust scattering and a pileup term (see text). The Suzaku 
and HETG spectra have not been renormalized with respect to one another. 
The residuals shown here are for this joint model fit with all broad and nan'ow 
line normalizations (both emission and absorption) set to zero. 

The above joint model received two modifications for 
solely the HETG s pectra. The simple_gpile2 model (see 
IHanke et al.ll2009l) was applied to describe the effects of pile 
up in the gratings spectra. Owing to this spectrum being 
fainter and harder than the Cyg X-1 spectrum described in 
IHanke et al.l ( 120091) . these spectra required an even smaller 
pileup correction. Additionally, we applied the dust scattering 
model (dust scat) to the HETG spectrum. As described in 
^ dust scattering represents a loss term for the high spatial 
resolution Chandra spectra. This component, albeit with a 
time delay of thousands to tens of thousands seconds, scatters 
back into the Suzaku spectrum from arcsecond to arcminute 
angular scales. To the extent that one can ignore the time evo- 
lution of the spectrum, the dust scattering term can be ignored 
in the Suzaku spectra. Lacking any detailed information of 
the average spectrum prior to the start of our observations, 
we do not apply any dust scattering correction to the Suzaku 
spectra. 

To be explicit, we apply to the Suzaku data a model of the 
form: 

TBnew* lines* 
(diskbb + powerlaw+ diskline + gaussian) (3) 

while we apply to the HETG data a model of the form: 

simple_gpile2 ® (TBnew * lines * dustscat* 
(diskbb + powerlaw+diskline + gaussian)) , (4) 

where lines represents the ionized line absorption. 
Throughout the rest of this work, we shall ignore as we are 
doing here the effects of dust scattering on the RXTE and 
Suzaku data. Again, this implicitly assumes that the bright 
phase spectrum has been steady over the time span of thou- 
sands to tens of thousands of seconds. This assumption is 
unlikely to be valid in detail; however, we lack the data to 
employ any more sophisticated assumptions. 
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Fig. 6. — The HEG data from a joint fit of tlie Suzaku and Chandra data 
with a simple continuum model (disk, powerlaw, broad and narrow Fe re- 
gion lines, neutral ISM absorption) plus two ionized absorber components 
described with the warmabs model. The residuals include all of these mod- 
els components; the remaining line residuals are features not described by the 
warmabs components. 

This simple model describes the spectra well with a dust 
halo optical depth of 0.25 ± 0.03 (90% confidence level), i.e., 
consistent with the fits to the color-color diagrams presented 
in S This spectral fit is shown in Fig. |5] The important 
point to note here is that the ionized absorption is extremely 
statistically significant in the Suzaku spectra, even though the 
individual lines are not resolved. Inclusion of these lines is vi- 
tal for obtaining a description of the soft end of the spectrum, 
and thus we include ionized absorption in all of our model fits. 
Being that we do not have HETG spectra simultaneous with 
our other three Suzaku spectra, we use the ionized absorption 
fit presented in Fig. |5] Specifically, for all subsequent spec- 
tral fits described in this work, we freeze the line positions, 
widths, and relative equivalent widths to the values found for 
this joint fit. The line equivalent widths are tied together via 
a single normalization constant, and this normalization con- 
stant becomes the sole fit parameter describing the ionized 
absorption. The model, however, still consists of 55 individ- 
ual absorption lines. As described below, this approach works 
well in all of our spectral fits. 

We explored whether a more sophisticated ionization model 
could be applied to our data. Rather than fit individual lines 
to the joint Suzak u and Chan dra data, we instead used the 
warmabs model jKallman & Bautista 2001) to describe the 
ionized absorption. A model that includes two ionized com- 
ponents — one with a column of (7 ±4) x 10^' cm"^ and ion- 
ization parameter logjo'C = 2.77 ±0.05 and one with a col- 
umn of (3 ± 2) X 10^' cm"^ and ionization parameter logj0 ^ = 
1.93 ±0.07 (90% confidence levels), and both with densi- 
ties of lO'^cm"-' — describes the data reasonably well. As 
shown in Fig. |6] however, a number of prominent line fea- 
tures are not well-fit in the Chandra data. Furthermore, aside 
from being an extremely computationally expensive model to 
run, when fit to the Suzaku data alone the warmabs mod- 
els would tend to gravitate toward low ionization parameter 
values that, although capable of mimicking some of the con- 
tinuum features of the Suzaku data, did not include the high 
ionization line features that were clearly present in the Chan- 
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Fig. 7. — The summed unfolded Suzaku spectra from unabsorbed periods 
from observations 2-A, shown in the Fe Kq line region. These data have been 
fit with a model consisting of an absorbed powerlaw, naiTow Fe Ka emis- 
sion line, naiTow Fe XXV absorption line, and a relativistically broadened 
diskline model. (See text for fitting procedures.) The individual model 
components (powerlaw, with added diskline. and separately with added nar- 
row emission and absorption lines) are shown overlain on the spectra, while 
the residuals are shown with and without the emission and absorption line 
components. 

dra data. For these reasons, we choose the empirical approach 
outlined above. 

Several things need to be borne in mind when considering 
our results: whereas we allow for an overall normalization 
change in our ionized absorber model, we do not allow for 
changes of the ionization state. As Fig. |5] shows, with the 
exception of a weak Fe XXV line, there are essentially no 
ionized lines with energies > 3 keV. Thus it is unlikely that 
this simple empirical approach will directly impact, for exam- 
ple, broad Fe line studies. However, in as much as ionization 
changes could affect our soft X-ray continuum fits, and these 
continuum fits in turn affect our estimate of the breadth and 
strength of the red wing of any broad line tail (see 35.31 ). po- 
tential ionization changes in this absorber must be considered 
an additional source of systematic uncertainty in the results 
discussed below. 

4.3. Composite Relativistic Line 

The fit described above consisted of both broad and narrow 
features in the Fe Ka region. Here we examine this spec- 
tral region in a somewhat phenomenological manner to gauge 
the interplay between these narrow and broad features. We 
wish to consider the average line profile from roughly simi- 
lar Cyg X-1 spectra. The three most recent Suzaku observa- 
tions occurred closest to orbital phase and show very similar 
color-color diagrams. As discussed below, these three spectra 
are fit with comparable ionized absorption normalization con- 
stants. We therefore consider these spectra jointly and ignore 
the first Suzaku observation. We fit these spectra with a model 
consisting of a powerlaw, narrow emission near 6.4 keV and 
narrow absorption near 6.7 keV, plus a relativistically broad- 
ened Fe line (again using the diskline model). The rela- 
tivistically broadened line is characterized by a broad tail ex- 
tending redward from the line rest frame energy, and a sharper 
peak (due to Doppler boosting) blueward of the line rest frame 
energy. 

For purposes of this discussion, we allow the powerlaw to 
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Fig. 8.— Suzaku (blue line), Chandia-HEG (grey line), and RXTE-PCA 
(purple diamonds) residuals from fitting an absorbed powerlaw to the 3.5- 
4.5 keV and 7-7.5 keV regions simultaneously. The Chandra-HEG data have 
been binned to match approximately the binning of the Suzaku data. 

assume independent parameter values for each of the Suzaku 
observations. The energies and widths of the narrow lines are 
tied together for all three observations, but the line strengths 
are fitted individually. The broad line has an energy fixed 
to 6.4 keV and an emissivity index fixed to /3 = -3 (i.e., line 
emissivity is c>c R~^, where R is the emission radius within 
the disk), and the disk inclination is fixed at 35° (the value 
adopted by Herrero et al.. 199 5, from the middle of the range 
of suggested incUnations). The inner radius of the broad line 
emission is tied together for all three observations, but the line 
strengths are left independent. The of the fit is determined 
from the sum of values from the individual observations; 
however, in Fig. |7] we show the spectra, fitted models, and 
residuals from the combined data. We fit the spectra only in 
the 4.5-7.5 keV bandpass. 

These joint fits strongly require the presence of a broad 
line. The broad line is significantly detected in each individ- 
ual XIS detector from each individual observation. Under the 
assumptions of the continuum model applied here, the broad 
line equivalent width is approximately three times that of the 
narrow line equivalent width. The best fit inner radius for 
the broad line emission is 10.6^; g GM/c^ (90% confidence 
level), i.e., close to the innermost stable circular orbit for a 
Schwarzschild black hole. As we will discuss further in ^5.11 
the specific value for such an inner radius is dependent upon 
the fitted continuum. Nevertheless, a broad excess is required 
redward of the Fe Ka line. 

In describing the residuals in this region with a broad line, 
we see that the Fe XXV absorption line occurs near the blue 
wing peak of the relativistically broadened line (Fig.|7]). Both 
the narrow emission and absorption lines in the Fe region are 
required in the HETG spect ra when considered by themselves 
(see also lHanke et al.ll2009l) . and thus they cannot be ignored 
in the Suzaku spectra. The HETG spectra, however, are not 
particularly well-suited for describing any broad component 
of the Fe line. (See Fig. |8] and the discussion of Hanke et al. 
l2009l concerning broad line fits to our prior joint RXTE- 
Ciandra-HETG observations of Cyg X-1.) 

On the other hand, if we compare the RXTE spectra to both 
the Suzaku and HETG spectra, we find that the RXTE spec- 




FlG. 9. — Possible X-ray emission geometry of the Cyg X-1 system. Part 
of the soft X-rays likely come from an accretion disk, which may be trun- 
cated at its inner edge by a hot corona. Two possible (of many suggested) 
geometries are shown for this corona: a torus, and a quasi-spherical cloud. In 
some models, the d isk can reform in the inner regions of the accretion flow 
IMaver & Pring le 2007). The corona will Compton upscatter a fraction of the 
disk component(s) soft X-rays to hard X-rays. Whether or not the corona en- 
velopes the disk(s) determines whether or not an additional, unComptonized 
disk component should be present in the spectra. (If the disk reforms on the 
inner edge of a toroidal corona, for example, a fraction of its soft X-rays 
will not intersect the corona.) Hard X-rays will be reprocessed into reflec- 
tion/fluorescent line components, predominantly by the outer disk. Finally, 
the jet might contribute to the observed X-ray emission via synchrotron and 
synchrotron self-Compton emission. 

tra alone require an Fe line equivalent width that cannot be 
accommodated with solely the narrow line component from 
the Suzaku-HETG fit. We highlight this fact in Fig. [8] where 
we show the residuals from a powerlaw fit in the 3.5^.5 keV 
and 7-7.5 keV region to the joint RXTE-Suzaku-HETG spec- 
tra. In this figure, we have binned the HETG spectra to match 
approximately the binning of the Suzaku spectra. There is 
overall good agreement between the RXTE and Suzaku resid- 
uals; both require a broad red wing in the Fe line region. The 
narrow Fe Ka emission and Fe XXV absorption are clearly 
seen in the HETG spectra; however, when viewed by itself, 
any broad Fe line component in the HETG spectra could be 
subsumed via a slight powerlaw slope change. We ascribe 
this latter fact to remaining calibration issues in the HETG 
detector. Based upon these results, we include broad and nar- 
row Fe Ka emission and narrow Fe XXV absorption in all the 
subsequent fits presented in this work. As we shall discuss 
below, however, the implied broad line parameters are depen- 
dent upon the assumed continuum model. 

5. BROAD BAND MODELS 

We now turn to a discussion of the 0.8-300keV, joint 
Suzaku-RXTE spectra for the bright phases of our four 
Cyg X-1 observations. Here we consider three different 
Comptonization models, two of which are discussed in detail, 
and a jet model. ("Toy geometries" for some of these situa- 
tions are presented in Fig. |9]) The latter model also contains 
Comptonization components that represent both Comptoniza- 
tion of disk photons as well as synchrotron self-Compton 
(SSC). The non-jet models discussed below all rely upon the 
eqpair model (Coppi 1999, l2005h to describe the fitted 
Comptonization components. 

5.1. Comptonization Descriptions 

The eqpair model allows for Comptonization by a coro- 
nal electron population that has both thermal and non-thermal 
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Fig. 10. — Unfolded Suzaku and RXTE spectra, plus model components and fit residuals from the high kTscd, thermal corona eqpair fits. Suzaku-XlS 
spectra from individual detectors/data modes and RXTE-HEXTE spectra from each cluster were summed for the figures, but not for the fits. Here, and throughout 
the remainder of this paper, XIS data are filled blue diamonds, PCA data are hollow purple circles, PIN data are hollow green diamonds, HEXTE data are brown 
squares, and GSO data are hollow magenta triangles. The following model components are shown: the summed model with the reflection component excluded, 
the reflection component by itself, the broad and nan'ow Fe Ka line components, the seed photons without Comptonization applied, the additional disk component 
with peak temperature tied to the seed photons, and the additional low temperature disk component. Neutral and ionized absorption are applied to all shown 
model components. 



energy distributions. The latter distribution is governed by a 
parameterized powerlaw, which for purposes of the fits de- 
scribed here we left at the default eqpair values. Specifi- 
cally, the electron phase space density follows a distribution 
oc E~^, where E is the electron energy, between energies of 
1 .3 eV-1 Me V. The eqpair model rather than being parameter- 
ized by a coronal electron scattering optical depth and tem- 
perature is instead parameterized by a seed electron optical 
depth (pair production within the corona, accounted for in the 
code, can yield higher net scattering optical depths) and coro- 
nal compactness parameters, £ (proportional to component lu- 
minosity or power divided by radius). The latter parameters 
are divided into two relative compactness parameters: £h/is, 
the coronal compactness divided by the seed photon compact- 
ness, and i„t/ih^ the compactness of the non-thermal electron 
population divided by the total electron population compact- 
ness. Over a wide range of parameters, the hardness of the 



spectrum and the temperature of the corona increases with in- 
creasing compactness (although pair production can modify 
that effect). 

The seed photon spectrum can be set to either a black- 
body or a disk spectrum (essentially the diskpn model 
from XSPEC; iGieriihski et al. 1999). We choose the latter 
seed photon spectrum, with the peak temperature of this disk 
(^Tdiski) becoming a fit parameter We also allow for an un- 
scattered disk component by including an additional diskbb 
component with peak temperature fixed to the peak temper- 
ature of the seed photons. Our choice for the normalization 
of this component depends upon the specific Comptonization 
model under consideration (see below). 

To account for soft excesses (i.e., below w 1.5 keV) in the 
spectra, we sometimes include an additional unscattered disk 
component with freely varying normalization and temperature 
(^disk2, kTdiski)- We Searched for, and easily found, solutions 
where the peak temperature of this unscattered disk com- 
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Fig. 1 1 . — Same as Fig. [To] but for the non-thermal corona, low kTsc^a eqpair model fits. The following model components are shown: the summed 
model with the reflection component excluded, the reflection component by itself, the broad and narrow Fe Ka line components, the seed photons without 
Comptonization applied, and the additional disk component with peak temperature tied to the seed photon peak temperature. Neutral and ionized absorption are 
applied to all shown model components. 



ponent was below w 300 eV. The eqpair model includes 
reflection from an ionized disk (i .e., a modification of the 
pexriv model; iDone et al]|1992h which is smeared by rel- 
ativistic distortions with the emissivity profile of a Shakura- 
Sunyaev type disk (Shakura & Sunvaev 1973). A fittable pa- 
rameter is the inner radius of this disk smearing profile (see 
below). We also fix the disk (i.e., reflector) temperature to 
10* K, the reflector inclination to 35°, and allow the reflection 
fraction (ft /2tt) and reflector ionization parameter to be fit 
parameters. The latter is limited to values < 1000. 

Line absorption is described by the parameteriz ed ( i.e., sin- 
gle normalization constant) model discussed in 34.21 A nar- 
row gaussian is added at 6.399 keV with its width frozen 
to 0.01 keV, but with a freely variable normalization. A 
diskline component (essentially t he profile for a d isk 
around a Schwarzschild black hole; IFabian et al.l 119891) is 
added to describe the broad line. Its energy is fixed to 6.4 keV, 
the disk line emissivity index is set to /3 = -3, and the disk in- 
clination is set to 35°. Typically one allows the inner emission 
radius to be variable (see 34.3l l. Here we tie this inner radius 



to that of the reflector in the eqpair model. The line nor- 
malization, however, is allowed to freely vary. 

A number of authors have considered "sphere-i-disk" Comp- 
tonization models wherein an inner quasi-spherical corona 
is encircled by a cool, geometrically thin disk that typically 
has a peak temperature of w 200 eV (iGierliriski et al.l 119971: 
iDove et al.iri997.) . We were unable to find any solutions for 
coronae with thermal electron distributions (i.e., int/ih frozen 
at 10"^) that allowed for such a low seed photon tempera- 
ture. The only such thermal corona solutions that we found 
required seed photon temperatures of 0.8-1 keV. Fit parame- 
ters for these solutions are presented in Table 4, and spectral 
fits are shown in Fig.fTOl 

Concepts for the potential geometry represented by this 
spectral fit are shown in Fig. |9] The high temperature seed 
photons could be a recondensed inner disk as envisioned by 
[Mayer & Pringle (2007), while the low temperature soft ex- 
cess, the broad line, and the reflection component could em- 
anate from an outer disk. (The narrow line likely arises from 
fluorescence from the secondary and/or spatially extended gas 
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Fig. 12. — Same as Fig.[To] but only showing observation 4 for the mul- 
tiple eqpair component fit (see text). The following model components 
are shown: the two individual eqpair components with their reflection com- 
ponents excluded, the reflection component of the dominant eqpair spec- 
trum, the broad and narrow Fe Ka line components, and the unscattered disk 
component. 

surrounding the system; Torrejon et al. 2010.) A more physi- 
cally self-consistent model would allow for an additional seed 
photon contribution from this outer disk, and possibly a broad 
line component from the inner, hot seed photon emission re- 
gion. However, it is interesting to note that these solutions re- 
quired very little additional unscattered, high temperature disk 
component. (Table 4 contains low- sometimes zero- values of 
^diski ■ Both the initial seed photon distribution and the associ- 
ated high temperature, unscattered disk spectrum are shown in 
Fig-ini the former is much more significant.) The amplitude 
of the soft seed photon contribution and the weakness of an 
additional unscattered high temperature disk would indicate 
a 'recondensed' disk region of only modest extent, predom- 
inantly within the confines of a more extended corona, from 
which we would expect little iron line contribution. 

On the other hand, these spectral fits do require an addi- 
tional, unscattered low temperature disk component. The val- 
ues of Aciisk2 presented in Table 4 correspond to inner disk radii 
of « 2-lQGM/c^, given a distance of 2.3 kpc, a black hole 
mass of 10 Mq, and an inclination of 35°. This is roughly 
consistent with the inner radius of the fitted broad Fe line 
and relativistically smeared reflector, for which we find val- 
ues ranging from Rin = 6-18 GM/c^. (The inner emission ra- 
dius of the line appears to be uncorrected with the inner ra- 
dius implied by the normalization of the low temperature disk 
component.) These radii values are roughly consistent with 
values close to the marginally stable orbit of a Schwarzschild 
black hole. 

The broad line amplitude is between 6-22 times the ampli- 
tude of the narrow line component. Fitted reflection fractions 
are f2/27r w 0.2. 

Comparable high se ed photon temperature fits were found 
by I Wilms et al.l 12006) when considering solely RXTE spec- 
tra. The major difference between those fits and the ones dis- 
cussed here is that our prior fits also required a significant con- 
tribution from an unscattered, high temperature disk. It is un- 
clear whether that difference is due to the inclusion of Suzaku 
data that extend to lower energies, the inclusion here of ion- 
ized absorption, or the fact that these spectra are for the most 



part fa inter and harder than those discussed by IWilms et alj 
( 120061) ■ The coronal compactnesses found here are compara- 
ble to the maximum values from our previous studies, while 
the seed optical depths are slightly lower (by « 0.3-0.5 for 
comparable observations). These values of seed optical depth 
and compactness correspond to coronal temperatures of 85- 
115keV and total optical depths of 1.1-1.4. The fitted reflec- 
tion fractions are slightly highe r here (by 0.05) than those 
fits with comparable ^h/^s fi'om lWilms et alj (120061) . 

We have found a set of Comptonization model solutions 
that do allow for a low seed photon temperature. For these 
fits, we include only one extra unscattered disk component 
(with its peak temperature tied to the peak seed photon tem- 
perature), and we further tie the inner disk radius to the inner 
radii of the eqpair reflector and the diskline emissivity. 
That is, we are explicitly mimicking the "sphere-ndisk" geom- 
etry (albeit not in a completely self-consistent manner) shown 
in Fig. |9] If we then relax the assumption of a purely thermal 
coronal electron distribution, good descriptions of the spec- 
tra are found with low seed photon temperatures. The ratio 
of non-thermal to total coronal compactness then falls in the 
range 4t/4 = 0.01-0.82. That is, up to 82% of the electron 
energy (presuming that the thermal and non-thermal electrons 
are cospatial) resides in a powerlaw distribution that extends 
to 1 MeV. Parameters for this model are presented in Table 5 
and spectra are shown in Fig.fTTI 

The quality of these fits is similar to or perhaps slightly 
better than that for the thermal corona fits. The fitted com- 
pactness parameters for the non-thermal coronae are approxi- 
mately d ouble those for the ther mal fits. (Similar results were 
found by llbragimov et alll2005l when comparing thermal and 
non-thermal eqpair fits.) At first glance, the optical depths 
seem smaller; however, owing to the large coronal compact- 
nesses pair production is significant and the net optical depths 
range from 0.7-1.6, while the thermal electrons have temper- 
atures that range from 55-160 keV. Reflection fractions are 
slightly larger (by up to 0.09) for these fits. The broad line in- 
ner radius in some cases has increased, and now ranges from 
w 6-32 GM/c^. Additionally, the amplitude of the broad line 
is now only 3-13 times that of the narrow line. This differ- 
ence compared to the thermal corona model is due to two fac- 
tors: the narrow line amplitude is slightly increased, while the 
broad line amplitude in two cases is nearly halved. 

We briefly consider one other thermal Comptonization so- 
lution that allows for a low seed photon temperature. (No "toy 
geometry" for this model i s shown in Fig . |9l) I t is a model 
similar to that considered by lFrontera et al. l ()2001|) when mod- 
eling BeppoSAX spectra of C yg X-1 and to that consid- 
ered bv lMakishima et al] (120081) when modeling Suzaku spec- 
tra of Cyg X-1. In these models two eqpair components 
are included. They share the same seed photon tempera- 
ture (tied to the peak temperature of an additional unscat- 
tered disk component, as before), and they share the same 
reflection parameters; however, they have independent com- 
pactness and optical depths. An example of such a model 
fit is shown in Fig. [121 These fits are generally successful, 
but no more so than the previous two models. Furthermore, 
owing to the larger degree of freedom given the two semi- 
independent coronal components, they are fraught with local 
minima representing qualitatively different relative contribu- 
tions of the two Comptonization components. We do not con- 
sider these more complex models further in this work, other 
than to point out that such double corona solutions as dis- 
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Fig . 13 . — The jet model fit to the estimated radio flux and measured X-ray 
spectra from observation 4. Individual model components are shown with 
both neutral and ionized absorption applied. We show the summed model, 
as well as the the jet synchrotron component, the jet SSC component, the 
Compton corona component, and the disk component. 

cussed bv lMakishima et alj (l2008h and iFrontera et alj (1200 lb 
can allow for both low seed photon temperature and purely 
thermal coronae. 

5.2. Jet Descriptions 

The Compton corona models described above provide a 
good description of the 0.8-300 keV spectra. These models, 
however, are not unique. Even within the restricted class of 
coronal models, we have shown three qualitatively different 
solutions that yield comparable fits. Furthermore, none of 
the Compton corona models give a self-consistent description 
of the correlated radio spectra. Although the four observa- 
tions discussed he re did not have simultaneous radio spectra , 
our prior studies doleissner et al] 120041: iNowak et aITl2005l: 
IWilms et al.ll200 6*) allow us to make good estimates of what 
the correlated 15 GHz radio flux likely was. We describe this 
estimated ioint r adio/X-ray spectrum with the jet model of 
iMarkoff. Now^& Wilms (2005J and Maitra et al.1 (IM09). 

It is important to distinguish the jet models as discussed 
in Markoff et al. (2004, 2005) and later works from the 
earlier iteration of the jet model discussed in Markoff et al. 
(2001, 2003). The 2001/2003 work considered a pai-ame- 
ter space where only synchrotron radiation from the jet con- 
tributed to the observed X-rays. Later work, and the version 
of the jet model used here, include X-rays produced not only 
by jet synchrotron radiation, but also by jet synchrotron self- 
Compton processes and Comptonization of seed photons from 
the disk. The synchrotron component is only mildly beamed, 
with peak /3 ~ 2 in the radio emitting portions of the jet. The 
inverse Comptonization component comes from near the jet 
base where (5 ^ 0.4, and thus undergoes very little beaming. 
Thus, although this does mean that the fitted jet parameters are 
dependent upon assumed inclination angle (e.g., the jet input 
power, A^j will depend upon assumed inclination), given the 
low /3 values inherent in the model it is possible to find good 
fit parameters for a wide variety of assumed inclinations. 

A full description of the main jet model parameters can be 
found in the Appendix of fMarkoff. Nowak & Wilmsl (l2005h 
and in iMaitra et al.1 (l2009h . The version of the jet model that 
we use is closest to that discussed in .Maitra et al.. (.20091) : 



therefore, we repeat part of the model parameter summary 
presented there. We assume a distance of 2.5 kpc and an in- 
clination of 35° for the Cyg X-1 system. The jet properties 
are then determined by the input jet power {Nj, expressed as a 
fraction of the Eddington luminosity for a IOM0 black hole), 
the electron temperature of the relativistic thermal plasma en- 
tering at the jet base (T^.), the ratio of magnetic to particle 
energy density (the equipartition factor, k), the physical di- 
mensions of the jet base (assumed to be cylindrical with ra- 
dius /?() and height /jq, the latter fixed to \.5R{)), and the lo- 
cation of the point on the jet (Zacc) beyond which a significant 
fraction of the leptons are accelerated to a powerlaw energy 
distribution that follows E~'\ with p being a fit parameter. 
The fit parameter eacc can be physically interpreted as being 
related to the particle acceleration rate (which is cx e^^'^,; see 
iMarkoff, Falcke & Fender 2001). 

As for the corona-only models, we also include emission 
from a low temperature disk with inner radius R\a and peak 
temperature Tdiski. (This disk is very similar in temperature 
and normalization to the low temperature disk component in- 
cluded in the Comptonization model fits.) A relativistically 
broadened Fe Ka line (with the same emissivity and inclina- 
tion parameters as for the corona-only models) is coupled to 
the disk component by having its inner radius of emission tied 
to the disk inner radius. A narrow gaussian line is included, as 
well as a reflection component (calculated with the reflect 
model from XSPEC). Unlike the reflection model internal to 
the eqpair code, the reflect model, which performs a 
convolution on any given input spectrum, does not account 
for relativistic smearing. The Fe edge near 7.1keV in this 
unsmeared reflected spectrum can produce a sharp feature not 
seen in the data residuals; therefore, we smear this component 
with a unit-normalized, cr = 1 keV gaussian convolutiorQ. 

Fit results for this model are presented in Table 6, Fig. [13] 
and Fig. [141 The overall quality of these fits is quite good, 
although not quite as good as for the corona-only models dis- 
cussed above. The jet models leave slightly larger residuals 
in the PC A spectra (near the Fe line region and the lOkeV 
region), as well as in the GSO spectra. The former could be 
related to inadequate modeling of the reflection spectrum (i.e., 
not using a relativistically smeared model). 

Fits to the hardest X-ray spectra (i.e., the cutoff seen in the 
GSO) are governed by the jet SSC and the inverse Comp- 
ton components. This latter component is comprised of a 
magnetized, beamed corona with a high temperature [« (3- 
5) X 10'" K, i.e., « 3-5 MeV], and low electron scattering op- 
tical depth (Tes ^ 0.01). A high temperature, low optical depth 
corona with a limited residence time under those conditions 
simplifies the jet code by allowing one to use a single scatter- 
ing approximation for calculating the Comptonization spec- 
trum. The compactness of such a corona, however, could lead 
to high pair production which would then serve to cool the 
corona an d increase its optical depth (see, for example, the 
critique of lMalzac. Belmont & Fabianll2009l) . 

^ Kernels for relativistic sme aring do ex ist, e.g., the models of 

IBrenneman & Reynolds! 420061) and lDauser et alj 420101) . The jet model cal- 
culations, however, are already substantially slower than the Compton corona 
calculations. Coupling them with the computationally expensive relativis- 
tic kernel would make the fits discussed here prohibitively time consum- 
ing to run. Furthermore, the jet model reflection geometry is likely more 
complex than that of the standard Compton corona: it should have separate 
contributions for the syn chrotron, SSC, and Comptonized disk photons (see 
IMarkoff & Nowakl2004l) . For these reasons we chose a very simple smearing 
profile. 
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Fig. 14. — Same as Fig. 1101 but for the jet model fits. The following model components are shown: the summed model with the reflection component excluded, 
the reflection component by itself, the broad and narrow Fe Kq line components, the jet synchrotron component (which dominates the 2-10 keV spectra for 
observations 1,2, and 4), the jet SSC component (which dominates the 100-300 keV spectra for observations 1,2, and 4), and the disk component. Neutral and 
ionized absorption are applied to all shown model components. 



Post-facto estimates of the p hysical self-co nsistency of fit- 
ted jet models are discussed by iMaitra et al.l (i200 9) who cal- 
culate the pair production and annihilation rates at the base of 
the corona. Fits to the hard state spectra of XTE J1181H-105 
yield far higher annihilation rates than production rates, and 
therefore represent a self-consistent solution. Fits to the hard 
state spectra of GX 339-4 yield mixed results. For some 
cases, annihilation and production rates are comparable and 
the resulting corona is marginally self-consistent. In other 
cases, the calculated production rates are an order of magni- 
tude higher than annihila tion rates, and th e resulting coronae 
are not self-consistent (M aitra et al.||2009 |). 

We have p erformed the same calculations as discussed by 
iMaitra et alj (l2009i) for jet models fitted in this work and find 
that all yield pair production rates larger than pair annihilation 
rates. The coronae within these jet models are therefore not 
self-consistent. Observation 3 produces a fit that is closest to 
being self-inconsistent, with pair production and annihilation 
rates being a factor sa 10 apart from one another. Jet fits to the 
other observations are further from self-consistency. 



Such inconsistencies were similarly true for early 
coronal models, e.g., t hose based upon the work of 
iHaardt & Maraschil (1199 Ih . A number of such "slab geom- 
etry" coronae, in order to produce the hardest X-ray spectra 
observed in sources such as Cyg X-1, used coronal temper- 
atures (apart from any consideration of pair production) that 
were unachie vable in those geometries (see th e discussion of 
this issue bv iDove. Wilms & Begelmanlll997l) . Such incon- 
sistencies were in fact what led to the development of more 
physically s elf-consis t ent co ronal models, e.g., the kotelp 
model of Stern et al. (1995) and iDove et alJ (Il997h and the 
eqpair model of Coppi (19991). 

We consider these jet models, however, as they are as of 
yet the only spectral models that make a serious attempt 
to explain the correlated radio spectra. Future iterations of 
these models will incorporate consideration of pair produc- 
tion, as well as synchrotron and SSC cooling of the Comp- 
ton corona (Pe'er & Markoff, in prep.). (As shown by 
iDove. Wilms & Begelman ,1997 and others, coronal cooling 
can be a more significant consideration than pair production 
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in some situations.) 

The derived parameters are similar to those found for fits 
to other BHC (see the fits to XTE Jl 1 18+105 and GX 339-4 
presented by Maitra et al. 2009). The power input at the base 
of the jet, in terms of Eddington luminosity, is comparable to 
the inferred bolometric luminosity (i.e., compare the values 
of A^j to the values of Table |2j although different inclination 
assumptions will alter the value of Nj, given the low /3 values 
in the jet, it is always comparable to the inferred system lu- 
minosity). The radius of the base of the jet is comparable to 
the radius of the innermost stable circular orbit (ISCO). The 
largest fit value exceeds the ISCO for a Schwarzschild black 
hole by only a factor of 1 1 . The shock acceleration distance 
along the jet, Zacc ranges from w 60-430 GM/c^. The former 
value is comparable to that found for XTE Jl 1 18+105, while 
the latter is comparable to that found for GX 339-4. As dis- 
cussed above, the fitted coronal electron temperatures at the 
base of the jet range from « (3-5) x 10"*K. Work in progress 
suggests that a physically self-consistent coronal temperature 
that properly accounts for the dominant effects of synchrotron 
and SSC cooling may truncate the allowed range of coronal 
temperatures to the lower end found here (Pe'er & Markoff, 
in prep.). Future versions of the jet code will explicitly ac- 
count for these effects. 

For all the fits presented in Fig. [14] both synchrotron and 
SSC radiation make substantial contributions to observed X- 
ray spectra. This is an important point to which we shall re- 
turn in ^ 

5.3. Comparison of Implied Relativistic Lines 

In Fig.[T5]we show the Fe line region residuals for the two 
Compton corona models and the jet model. These residuals 
are derived by setting the line normalizations (including ab- 
sorption line normalizations) to zero, without refitting. They 
are also different from line residuals that are often used for 
illustrative purposes where only a portion of the continuum 
outside of the line region is fit with a simple powerlaw (e.g.. 
Fig. IT). The residuals in Fig. [15] are meant to highlight the 
differences in the derived line parameters for the three main 
continuum models discussed here. 

As shown by the parameters in Tables 4-6, there are in- 
deed differences among the inferred line parameters. For the 
most part, the purely thermal Comptonization model yields 
the strongest line equivalent widths, while the non-thermal 
Comptonization model yields the weakest equivalent widths 
(by up to a factor of two). The thermal Comptonization model 
usually also yields the smallest inner radii for the relativisti- 
cally broadened line, in some cases extending all the way to 
the ISCO for a Schwarzschild black hole. Larger values of 
these radii, however, are typically found with the other two 
models. These differences are evident in the lower energy 
residuals shown in Fig. [7] 

To be clear, a substantially broadened line is required by 
all spectral models that we have explored, and there is a good 
agreement between the PC A and Suzaku-XlS residuals re- 
garding this fact. We also reiterate that for all explored mod- 
els, the equivalent width of the fitted narrow line (typically 
< 15 eV) is only a fraction of that for the broad line, and is in 
fact substantially smaller than the Fe line equivalent width we 
have found for any mod el that we have app lied to any RXTE 
observation of Cyg X-1 dWilms et alj|2006f) . Furthermore, for 
these four observations, we have not found any implied inner 
radius greater than w 40GM/c^. The precise details of this 



broadened line, however, do depend upon the assumed con- 
tinuum model. 

6. DISCUSSION 

In this work we present broad band (0.8-300 keV) fits to 
four separate observations of Cyg X-1 that have simultaneous 
Suzaku and RXTE data. The most recent of these observa- 
tions also has simultaneous Chandra-HETG data. Each of 
these observations shows evidence of dipping events likely 
due to dense structures ("clumps") in the otherwise highly 
ionized wind of the secondary. This is seen in the lightcurves 
(Figs. [TJ, and the color-color diagrams (Fig. [2] which can 
be modeled by the presence of a dust scattering halo). The 
ionized absorption is very significant in the Suzaku spectra 
(Fig. O even though in this work we consider the least ab- 
sorbed periods of the lightcurve. 

• Modeling the ionized line absorption present in Cyg X- 
1, via the use of parameterized fits to the HETG spectra, 
is crucial for deriving good fits to the soft X-ray spectra 
seen by Suzaku and RXTE. 

We note that even accounting for this ionized absorption, 
the spectral fits presented in this paper yield reduced that 
range from 1.8-2.5. Do such values truly represent good fits 
to these data? Cyg X-1 is bright, and these observations are of 
sufficient length that the signal-to-noise values for these spec- 
tra are quite high. The spectra are dominated by systematic 
errors, especially at the soft X-ray energies. We already have 
added 0.5% systematic errors to the PCA spectra, which is a 
reasonable estimate for the internal uncertainty of the PCA. 
There are also relative uncertainties among the detectors. For 
observation 4, if we increase the PCA systematic errors to 1 %, 
add 1% systematic errors to the HEXTE A cluster (i.e., the 
fixed cluster), and add 3% systematic errors to the Suzaku- 
XIS spectra, then the fits presented here would have reduced 
sa 1 . Comparing the fits among the individual XIS spectra 
(representing both different individual detectors and different 
data acquisition modes), we have found that ±3% is a rea- 
sonable estimate of the end-to-end differences among these 
spectra. We hypothesize that the quality of the fits presented 
here is near "optimal" given the current internal and relative 
calibrations of these detectors. 

The observations that occurred on 2007 May 17, i.e., the 
third set of observations, are potentially the most problem- 
atic in turns of cross-calibration issues. For our other sets 
of observations, the RXTE-PCA data act as a "bridge" be- 
tween the Suzaku-XlS and -PIN spectra. Each can have its 
relative normalization anchored by a comparison to the PCA, 
which overlaps the energy coverage of both detectors. This 
is lacking in the third observation, leading to the worry that 
changes in the normalization constant are subsuming, for ex- 
ample, continuum spectra associated with the spectral break 
at « lOkeV. For the third observation, we find the ratio of the 
PIN to XIS 1 normalization constant to range from 1 .10-1 .17. 
This is to be compared to the 1.16-1.18 value found for fits to 
the first observation, the 1.06-1.08 found for the second ob- 
servation, and the 1 .10-1 . 1 1 found for the fourth observation. 
The expected valu^l for the PIN/XIS comparison is 1.16. 
There is some amount of scatter for the cross-normalization 
values among the different observations; however, the first, 

^| http : //hea sarc ■ gsf c ■ nasa ■ gov/docs/ suzaku/ anal ysis?] 
watchout . html . Since we did not fit XlS data for the fourth observa- 
tion, here we compare to XIS 1 . 
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Fig. 15. — Fe Ka line region residuals for observation 4, obtained from our best fits by setting the broad and narrow line normalizations (including the 
normalization of the Fe XXV absorption line) to zero. Left: the high ^cTsecd eqpair fit, the hybrid thermal/non-thermal eqpair fit, and the jet model fit. (Purple 
circles are the RXTE residuals, while the blue histograms are the Suzaku residuals.) 



second, and fourth observations show little scatter in its value 
for different fits to the same data. More scatter is seen for 
the third observation; therefore, additional systematic uncer- 
tainties need to be considered as being present for that set of 
observations. 

Using simple broken powerlaw and exponentially cutoff 
powerlaw fits, we find that these spectra are among the hardest 
seen in the "low hard state" of Cyg X-1 over the past decade. 
For all four observations, the spectra are clearly detected out 
to 300 keV with the GSO, and exponential cutoffs are well 
constrained. 

• Although these are among the hardest Cyg X-1 spectra 
ever detected, the exponential folding energies vary by 
over a factor 1.5, and range from 160-250keV. 

Historically observed hard state spectra in Cyg X-1 have 
shown folding energies that vary over a slightly wider range, 
while hard state BHC as a class show folding energies that 
span a factor of 5. 

As these spectra are among the faintest and hardest for 
Cyg X-1, they make excellent test beds for theoretical models 
that posit, for instance, that the hard state represents a con- 
figuration with an inner disk that has evaporated into a quasi- 
spherical corona. In such a scenario, we might expect these 
spectra to show the most "extreme" recession of the inner 
disk, although we have noted that the bolometric luminosi- 
ties represented by these spectra only span a factor of two. At 
a few percent of the Eddington luminosity, they are not far 
below the expected soft-to-hard state transition. Numerous 
transient BHC sources show much fainter hard states as they 
fade into quiescence. 

We have presented a number of different spectral models, 
all of which describe the 0.8-300 keV spectra well. Some of 
these models describe the X-ray spectra primarily with Comp- 
tonization components (whether due to a thermal or hybrid 
thermal/non-thermal corona), while the jet model is domi- 
nated by synchrotron and SSC emission from the jet. All of 
these models have a number of features in common. 

• All models require a soft excess that here we describe 
with a disk component with low (A:rdisk ~ 200 e V) peak 
temperature. The implied inner radii of these disks 
range from 2-40 GM/c^. 

• All models require a relativistically broadened line 
component. The inner emission radius of this broad- 
ened line never exceeds w 40 GM/c^, but for some 
models is as low as 6 GM/c^. 



• All models require a reflection component. The typical 
values for the reflection fraction are H/Itt w 0.2-0.3. 

• All models imply that the spectral hardening at w 
10 keV is not solely due to reflection. 

This latter point is very important, and broadly agrees 
with similar co nclusions dr awn b y Frontera et al. (200f), 
llbragimov eTaTI 62005 ) . and iMakishima et al.. (2008). The 
presence of the broad Fe Une and some of the spectral curva- 
ture in the 20-300 keV band is a clear indication of the pres- 
ence of reflection. However, this reflection spectrum is not sit- 
ting on top of a simple "disk+powerlaw" spectrum. There 
is additional continuum complexity separate from reflection 
that contributes to this perceived break. A high seed photon 
temperature in the thermal corona model yields a soft excess 
in the 2-lOkeV band and thus contributes to the measured 
break in that scenario. (See the discussion in Wilms et al] 
I2OO6.) As has been discussed by Ibragimov et al. (2005), a 
non-thermal electron population in the corona can lead to a 
soft excess in the 2-lOkeV band that helps contribute to what 
otherwise would be modeled as a r eflection break at lOkeV . 
The two corona model discussed bv'Ma kishima et all (120081) . 
a version of which is shown in Fig. [121 rather explicitly re- 
places part of the reflection component with a broad band 
continuum model. Finally, in the jet paradigm the spectral 
break at lOkeV is partly attributable to the transition from 
dominance by synchrotron emission to SSC emission in the 
continuum. 

There are plausible physical scenarios for each of the 
discussed spectral models, with some hypothesized geome- 
tries being shown in Fig. |9] The hybrid thermal/non- 
thermal coronal model is the closest to the concept 
of the quasi-spherical inne r corona with outer geo- 
metrically thin disk (e.g., 'Eardlev, Lightman & Shaiprd 
1975; Shapiro. L ightman & Ea rdlev 1976; Ichimaru 19771 
iDove et al.l il997. etc.). We have only been able to find such 
solutions, however, when invoking a partly non-thermal elec- 
tron population in the corona. 

The purely thermal coronal model contains two disk com- 
po nents that are reminisce nt of the physical description given 
by iMayer & Pringl^ (2001). These authors describe a situa- 
tion where an outer, geometrically thin, cool disk surrounds 
an inner, geometrically thick, hot corona. In the very inner 
radii of this corona, however, thermal conduction leads to it 
condensing into a geometrically thin and optically thick disk 
(see Fig. |9|l. Such a component could supply the high tem- 
perature seed photons in our thermal corona solutions, while 
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the lower temperature outer disk could provide the bulk of the 
reflection features. The jet model has a natural physical in- 
terpretation in that the usually observed optically thick radio 
spectrum observed in the hard state is clear indication of the 
presence of a jet. The question that remains is the contribution 
of this component to the X-ray band. 

The fact that the continuum is more complex than a simple 
"disk+powerlaw", yet there are multiple, physically mo- 
tivated models that yield comparably good spectral fits, leads 
to the final point. 

• Although a relativistically broadened line is required in 
all of our spectral models, the parameters of this line 
are dependent upon the presumed continuum model. 

Coupled with this dependence upon assumed continuum spec- 
trum is an implicit dependence upon ionized absorption, for 
which we have detailed Chandra-HETG measurements for 
only the fourth observatioifl Again, we have not found an 
inner radius for this line that exceeds « 40 GM/c^. Given the 
variations of this line with presumed continuum, however, we 
are as of yet unable to use this line for more refined diagnos- 
tics such as estimates of black hole spin. 

Although we are unable as of yet to draw firm conclusions 
as to the best geometry and physical mechanisms to describe 
the hard state spectra of Cyg X-1, these new joint Suzaku- 
RXTE data provide a st unning contrast to our prior results 
using solely RXTE data (IWilrns et al.ll2006^ . For the RXTE 
data alone we were able to describe the 3-125 keV spectra 



with a variety of physically motivated Comptonization mod- 
els and to describe correlations among the fit parameters. On 
the other hand, the simple exponentially cutoff, broken pow- 
erlaw models with a single, bro ad gaussian Une d escribed the 
data equally well, if not better (Wilms et al. 2006). When con- 
sidering the 0. 8-300 keV Suzaku-RXTE data discussed here, 
this is no longer the case. We now require complex absorption 
at low energy, an asymmetric broad line plus narrow emis- 
sion and absorption components in the Fe line region, and a 
complex continuum model. The catalog of Suzaku observa- 
tions of BHC will continue to increase such that we observe 
a wider variety of BHC states and luminosities. Furthermore 
the sophistication and physical self-consistency of the spec- 
tral models will continue to improve. Together, they offer the 
promise of obtaining a better understanding of the physical 
processes occurring in these BHC systems. 
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APPENDIX 

Suzaku ATTITUDE CORRECTION AND PILE-UP ESTIMATION 

We have created two tools to aid our Suzaku data analysis: aeattcor . si and pile_estimate . si. The tools, descrip- 
tions of their use, and example results can be found at: http : //space .mit . edu/ASC/software/suzaku[ The former 
tool further corrects the Suzaku attitude solution, whereas the latter estimates the degree of pileup in a given subset of the obser- 
vation. Both are scripts written in^S^Lang but are designed to be run on a terminal command line using ISIS, the Interactive 
Spectral Analysis System (iHouck & Denicol a 2000), as a driver. For all intents and purposes, the scripts behave similar to the 
typical Unix command line tools found in the HEASOFT package. 

Thermal flexing of the Suzaku spacecraft leads to a slow wobbling of the optical axis, and hence blurring of the image. 
Current Suzaku tools partiall y correct this ef fect by adjusting the spacecraft attitude file based upon details of the spacecraft 
orbit, temperature, etc. (Uchiv ama et al.ll2008l) . aeattcor . si further improves this correction (Fig.[T6]l for a bright source by 
using its time-dependent detected image to create a new attitude file. Specifically, the tool attempts to shift the time-dependent 
mean detector image position to a new specified, fixed sky position. It presumes that there are no intrinsic variations in the time- 
dependent image position. (A mean image position may be time varying, for example, if it is comprised of two or more variable 
sources of comparable flux. Additionally, a highly variable and piled up source where an "image crater" appears and disappears 
over the course of the observation may lead to a variable mean image position.) The bin time over which the tool searches for 
and shifts the image peak is a user selectable parameter; however, we have used the default of 100 seconds. 

We applied the tool to detector images that first underwent the standard attitude correction process (UchiyamaetaDHOQl). We 
chose circular regions with radius « 2' that were visually centered on the source image. The tool then created a new attitude 
correction file, which we applied with the xiscoord tool. All of our attitude correction files were created using solely the XIS 
1 images. 

Pileup occurs in CCD detectors when two or more photons fall on the same or neighboring pixels dur ing th e same 
readout frame, and therefore are read as a single higher energy photon or are discarded as a bad event jDavisI |2001). 
pile_est imate . si was run after aeattcor . si, as uncorrected blurring leads to an underestimate of the degree of pileup. 

The tool first creates a lightcurve in three energy bands, using both rate and counts. Here we used the default values of 0.5- 
l.SkeV, 1.5- 3 keV, and 3-9 keV for the energy bands and 32 seconds for the time bins. These data are then passed to the vwhere 
(lNoblell2005h filtering tool. The vwhere tool was used to create the color-intensity diagrams shown in Fig.|2]i, and it was used 
to select times of approximately uniform rates and colors (see These time selections were then written to a filter file used for 
subsequent data extraction using the xselect tool. 

A box-car smoothed (3x3 pixel bins) image is displayed using the ds 9 tool, with the image rescaled to the estimated pile up 
fraction. Pile up fraction here is defined to be the ratio of events lost via grade or energy migration to the events expected in the 
absence of pileup. Furthermore, the pile up fraction is based upon the mean counts per 3x3 pixel region per readout frame. The 
image displays discrete steps that represent the minimum pileup fraction in the displayed region, with the exception of the highest 
displayed value. This value corresponds to the image maximum, and is shown over regions at this maximum down to half way 
towards the next lowest displayed value. 

Using this image, we chose to exclude the central regions with pile up fractions 10% (Fig.fTSI). Typically, this was an « 20" 
radius region, comprising w 1/3 of the detected events. The pi le_est imate . si tool then gives an estimate of the mean level 
of the remaining pile up fraction in the image, which for these data was typically < 4%. 
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Fig. 16. — Left: An example of Suzaku attitude coiTection. The left half of the image shows a Cyg X-1 image using the standard attitude correction. The right 
half of the image shows the improvement with aeattcorr . si. Right: An image of Cyg X-1 where discrete colors correspond to the pileup in that region. The 
outer white circle denotes the outer boundary of an annular extraction region, and the inner white circle denotes the inner boundary. Pileup fractions within this 
excluded region are as high as 35%. The average effective residual pileup level is < 4%. 
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Parameters for thermalComptonization spectral fits 
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Note. — See text for a description of the model. Spectra are normaized to PCA data. Disk normalizations are ^in/D)^ cos 9, where Ria is the inner disk radius in units of km, D is the source distane in units 
of lOlqK, and is the inclination. Line normalizations are integrated phton flux in the lines (y/cm^/s). Error bars are 90% confidence level for one interesting parameter. 



TABLE 5 

Parameters for non-thermal Comptonization spectral fits 



Date 
(yyyy-mm-dd) 


Nh 

(10^2 ^-2^ 






4t/4 




n/2-n- 


i 

(47rfi„„/«) 


(GM/c^) 


(keV) 


''^skline 

(10-3) 


^line 

(lo-*) 


Glines 


X^/DoF 


2006- 10-30 

2007- 04-30 

2007- 05-17 

2008- 04-19 


081!S:g| 

86+° '" 
■ -Q.Q2 

'■OO-O.OI 


7 77Q+0.027 
^■^^'-0.082 
3 158+O.003 

J 811+0.383 

1 1^R+^:8I? 


16.04!°;°^ 
13.29!°! 
12.98:*:| 
15.28^:g 


n C7+0.0I 

A 01+0.01 
"■"'-0.01 

0.36!°;°' 
"■783:g 


1 ^9+0.03 
'■■"-0.02 
40+0.01 
^■^^-0.06 

01+° * 

1 4i3i.0I 

-0.02 


21 1+0.001 
"■^"-0.008 

305+0.010 
0257+0038 
239+'^® 

"■■^-^^-0.002 


499!^J 
97+f 

0+18 

189!^? 


6 7+0.1 

"■'-0.1 

6 3+20 

'I'l t+LS 


n 1Q9+0.002 
■J-A^^-O.OOl 
0.158!°;°°! 

0.150^* 


"■"-0.9 
3 0+0.6 
^■"-0.6 

3 4+0.6 

4 1+^:^ 

^■'-1.0 


1 1 1+3' 

"■'-2.9 

9.7;[-' 
2.6+' 8 
9 4^:1 


32+° °'' 
"■^^-0.06 

n ss+0.02 

U.S5_„„2 

81+* »' 

"■0^-0.09 
'•""-0.07 


lliABll 
1919/844 
810/485 
1535/845 



Note. — See text for a description of the model. Spectra are normaized to PCA data. Disk normalizations are ^in/O)^ cos 0, where Rya is the inner disk radius in units of km, D is the source distane in units 
of lOlqK, and Q is the inclination. Line normalizations are integrated phton flux in the lines (cf/cTc^li). Error bars are 90% confidence level for one interesting parameter. 



TABLE 6 

Parameters for jet-dominated spectral fits 



Date 


Nh 




«0 


Zacc 


Cacc 


P 


k 


7; 




^diskl 


n/2v 


-'^skline 


^line 


Glines 


X^/DoF 


(yyyy-mm-dd) 


(10=2 cm-2) 


(10-2 tEdd) 


(GM/<^) 


(GM/c^) 


(10-*) 






(10'°K) 


(GM/c^) 


(10*' K) 




(10-') 


(10^) 






2006- 10-30 

2007- 04-30 

2007- 05-17 

2008- 04-19 


0.86!°;°! 
1 15+0.01 

96+'^:'^' 

"•'"-0.01 

1 00+°°' 

'•""-0 01 


3 07+" "' 
-••^'-0.07 

3 22+0.05 
1 02+8:81 

' ^"^-0.05 
o o^+O.Ol 
^•"■'-0.02 


9 90+0.01 

^•^°-o.oi 

9 97+0.01 
• -Q*OJ 
77 2+2 

2.00+^^ 


66.9:j',' 
428!' 1° 


89+° °' 
"•"'-0.01 

1 43+0.01 

1 62^^^ 

'•"■'-0.02 

73+<'0i 
"•'^-0.01 


1 ns+o.oi 
^■"8-0.01 
1 79+0.01 

9 33+8:81 

■^ -0.01 

1 70+°°' 

'■'"-0.01 


1 75+0.30 
'■'^-0.15 

in+° °i 

"■'"-Q.Ol 
4 54+0.16 
^■■^^-0.09 

39+0 01 

"■^'-001 


-> 97+0.03 
^■^'-0.01 
5 15+0.01 

4 14+8:81 

^■'*-0.02 
71: +0.05 
^- -0,06 


7 1+°' 
'•'-0.5 

28.4+»-* 
15 5+'':^ 

99 Q+0.4 
^^■'-0.6 


2 34+0.03 
^■-'^-0.01 

1 66+°°2 
'■""-Q.Ol 

1 79+0.01 

'■'^-0.01 

■■86!|j:Sl 


1 1+°°' 
"•"-0.01 

3()+° °2 

34+°-™ 
"•^^-0.01 

24+° °' 
"•^^-0.01 


1 1 2+°-' 
"•^-0.8 
3 9+0.5 

'••^-0.5 

6 9+°-* 

"•'-0.9 


1+2,0 
7 7+2.0 

34rf:9 
^•^-1.5 

8.9!*;] 


A -50+0. 06 
"•■'''-0.03 
0.84!°;» 
0.46!«;i 


2444/976 
2042/843 
831/484 
1543/844 



Note. — See text for a description of the model. Spectra are normaized to PCA data. Line normalizations arc integrated photonflux in the lines (7/cm^/s). Error bars are 90% conftdcncc level for one 
interesting parameter. 



